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Abstract 

Continuum models with critical end points are considered whose Hamiltonian H[<f), ip\ 
depends on two densities (J) an d ip. Field-theoretic methods are used to show the 
equivalence of the critical behavior on the critical line and at the critical end point 
and to give a systematic derivation of critical-end-point singularities like the ther- 
mal singularity ~ |i| 2_Q of the spectator-phase boundary and the coexistence sin- 
gularities ~ or ~ \t\@ of the secondary density (ip). The appearance of a 
discontinuity eigenexponent associated with the critical end point is confirmed, and 
the mechanism by which it arises in field theory is clarified. 
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1 Introduction 



Critical end points are ubiquitous in nature. They occur when a line of criti- 
cal temperatures T c (g), depending on a nonordering field g such as chemical 
potential or pressure, terminates at a line ga{T) of discontinuous phase tran- 
sitions [1-3]. Two familiar examples are: (i) the critical end point (CEP) of a 
binary fluid mixture where the critical line of demixing ends on the liquid-gas 
coexistence curve; (ii) the CEP of 4 He, the terminus of the lambda line on 
the gas-phase boundary. On the critical (or lambda) line the disordered and 
ordered phases separated by it become identical critical phases; in the case of 
a binary fluid with components A and B, the disordered phase corresponds to 
a homogeneously mixed fluid af3, and the ordered ones to an A-rich phase a 
and a B-rich phase /3; in the case of 4 He, the disordered and ordered phases 
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are normalfTuid and superfluid, respectively. A crucial feature of a CEP is that 
a critical phase coexists with a noncritical ('spectator') phase 7 there. 

Although CEPs were encountered in numerous studies of bulk and interfacial 
critical phenomena [1-11, and their references] in the past decades, they have 
rarely been investigated for their own sake. This may be due to the expectation 
that the critical phenomena at a CEP should not differ in any significant way 
from critical phenomena along the critical line T c (g) [1]. However, recently 
it has been pointed out [2-4] that even the bulk thermodynamics of a CEP 
should exhibit new critical singularities, not observable on the critical line. On 
the basis of the phenomenological theory of scaling it was predicted [2] that 
the first-order phase boundary g a (T) should vary near the CEP (T = T e ,g = 
9c = ga(T c )) as 

9<t(T) - g^{T) - -jS |t| 2 "« (1) 
(2 — a)(l — a) 

in the limit t = (T — T c )/T c — > 0±, where (y£ eg (T) is regular in T. Furthermore, 
the amplitude ratio X+/ X°_ should be equal to the usual universal (and hence 
g independent) ratio A + /A_ of specific heat amplitudes A±. These are defined 
by writing the specific heat singularity at constant g 7^ g e on the critical line 
as A±(g) \T — T c (g)\~ a . In other words, the singularities displayed by g a {T) 
should be of the same form as those of the bulk free energy of the disordered 
(a/3) and ordered (a + j3) phases near T c (g). 

The phenomenological scaling arguments leading to (1) can be extended in a 
straightforward fashion to determine the singularities p g , the thermodynamic 
density conjugate to the nonordering field g, should display as the CEP is 
approached along the coexistence boundary [12]. They yield the singular part 

pf s = p 9 - p7 s ( t ) « u °± \tf + vg \tt a . (2) 

Having in mind binary fluid mixtures, we take the CEP to lie on the liquid 
(rather than the gas) side of the coexistence boundary. The quantity p g may be 
identified as the total density of the fluid. For a hypothetical symmetric binary 
fluid whose properties are invariant with regard to simultaneous interchange 
of its two constituents A and B and their respective chemical potentials /^a 
and hb, the amplitudes U± would vanish. More generally, this would be true 
for systems that are describable by a continuum Hamiltonian which is even 
in the order parameter field 0. Just as X^/X ^, the ratios as well as 

U+/U°_ (if C/° 7^ 0) are universal and can be expressed in terms of standard 
universal amplitude combinations [12]. 

The |t| 2 ~ a singularity of (1) has been checked by Monte Carlo calculations [12] 
and verified for exactly solvable spherical models [8,9]; the \t\ l ~ a singularity 
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of (2) is consistent with the jump in the slope of p g {T) found in mean field 
and density functional calculations [10,11] and has also been seen in Monte 
Carlo simulations [12,13]. 

Here we will address the issue of CEP singularities via the field-theoretic renor- 
malization group (RG) method. This approach is known to provide both a con- 
ceptually reliable basis of the modern theory of critical phenomena as well as 
powerful calculational tools (see, e.g., [14,15]). Surprisingly, it has not yet been 
applied with much success to the study of CEPs. We are aware of only one such 
work that goes beyond the Landau approximation, an (e = 4 — d)-expansion 
study of a scalar 8 model with negative 6 term [16]. Its one- loop result is 
that the critical line and the CEP are controlled by the same, standard 0(e) 
fixed point. Unfortunately, the model investigated has rather special features: 
its first-order line does not extend into the disordered phase; as its CEP is 
approached from the disordered phase, the order parameter (0) becomes crit- 
ical and exhibits a jump to a nonvanishing value upon entering the ordered 
phase; and no critical fluctuations occur in its ordered phase. Hence it clearly 
does not reflect the typical CEP situation in which the two-phase coexistence 
surface bounded by the critical line T c (g) meets the spectator phase boundary 
in a triple line; its applicability appears to be quite limited. 

One should also note that the above RG scenario differs from the one found 
in position-space RG calculations [5,6] of lattice models with conventional 
CEPs. In the latter scenario the critical line and the CEP are mapped onto 
separate fixed points, where the CEP fixed point has two relevant RG eigen- 
exponents that are identical to those of the former, plus the additional one 
y = d, characteristic of discontinuity fixed points [17], but absent in [16]. 

We conclude that systematic field-theoretic RG studies of appropriate models 
are urgently needed. A first obvious goal one would hope to achieve is a sys- 
tematic derivation of the singularities in (1) and (2). This involves showing the 
equivalence of critical behavior at the CEP and on the critical line .[] Provided 
the above RG scenario with two separate fixed points prevails, one must prove 
that the associated critical spectra match, demonstrate the existence of the 
discontinuity eigenexponent y = d and clarify its significance. 

We have recently carried out such an investigation. In the sequel, we will briefly 
describe the main steps of our procedure and our findings. A more detailed 
exposition of our work will be given elsewhere [18]. 



1 In their excellent survey [3] of the present state of the theory of CEP singularities, 
Fisher and Barbosa warn that this equivalence, with matching critical spectra of the 
corresponding two fixed points, need not be an invariable rule, even though their 
work confirms it, just as our own. 
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2 Models 



First, we must choose an appropriate continuum model. Natural candidates 
are models whose Hamiltonian H[(f),ip] depends on two fluctuating densities: 
a (primary) order parameter field <f>(x) and a secondary (noncritical) density 
ip(x). The form of Ti can be guessed on purely phenomenological grounds, 
but can also be derived by starting from an appropriate lattice model, such as 
the Blume-Emery-Griffiths model [19] on a rf-dimensional simple cubic lattice. 
This is a classical spin S — 1 model with Hamiltonian 



n BEG [s] = - £ [J Si Sj + k si s* + l (si Sj + s t sf) 

(i,3) 

-J2{HSi + DSl), Si = 0,±1, 



(3) 



where indicates summation over nearest-neighbor pairs of sites. We pre- 
sume the interaction constants K and J to be positive ('ferromagnetic'), and 
L > 0. The quantities H and D correspond respectively to even and odd linear 
combinations of the chemical potentials /xa and /ib [20]. 

Performing a Gaussian ('Kac-Hubbard-Stratonovich') transformation with re- 
spect to both {Si} and {Sj}, one can map the model (3) exactly on a lattice 
field theory with fields (pi G R and ipi G R [18]. To make a continuum approxi- 
mation, we replace these by smoothly interpolating fields 4>{x) and ip(x), and 
Taylor expand nearby differences (pi — (f)j about their midpoint (i + j)/2. We 
thus arrive at a continuum model with the Hamiltonian 



Hi[<t>\-- ld d x 



n 2 



H 



12 



d d x 



2 



ofx 



do 



(4) 
(5) 

(6) 
(7) 



The 3 and ip s terms have been eliminated by shifts (p(x) — > </>(cc) + O and 
V'(aj) — ► ip(x)+tp . Monomials of higher order and higher-order gradient terms 
have been dropped. 

The terms retained in 7i require explanation. Consider, first, the case of a 
symmetric CEP, in which H[—(p,ip] = Tl[(p,ip], i.e., h — du — en = 0. Owing 
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to this symmetry, <f> and if) do not 'mix' and hence may be chosen as the fields 
that become critical or remain noncritical at the CEP, respectively. To assess 
the relevance of contributions to 7i via power counting, the coefficients A and 
b 2 should be taken dimensionless. Thus <fi has the usual momentum dimension 
[(f)] = (d — 2)/2, while [if)] = d/2. In [</>], we have kept all monomials of the 
standard 4 Hamiltonian, namely those (except <f> 3 ) having coefficients with 
nonnegative momentum dimensions for e > 0. For the remaining interaction 
constants, one finds [g] = — [e 2 i] = 2 — e/2, [d 2 i] = e/2, [B] = —2, and 
[64] = e — 4. This suggests that B, e 2 i, and b 4 may be expected to be irrelevant 
in the RG sense and hence can be set to zero. If we did this, 7i would reduce 
to the Hamiltonian of the dynamic model C [21]; it would be quadratic in 
if), so if) could be integrated out exactly. The resulting effective Hamiltonian 
would be identical to Hi [<f>], up to a change of its parameters a 2 and a 4 , and 
an overall constant. 

The terms oc B, en, and e 2 i have been introduced because they play a role in 
the analysis of inhomogeneous states with a liquid-gas interface^ [18]. Since 
our main focus here is on bulk critical behavior, we can indeed set B = e 2 i = 
in the sequel. However, 64 must not be set to zero because, then, we would 
not be able to describe a/5-7 coexistence, nor would the model have a CEP. 



3 Landau theory and beyond 

Application of the Landau approximation to the model (4)— (7) yields a phase 
diagram with a CEP and the correct topology, provided its parameter values 
are in the appropriate range [10,13,18]. With the choices a 2 < 0, d 2 i > (aside 
from a 4 > and 64 > 0), one finds a critical line with 

rl 2 h n 2 

W = 4>c = -02/021 , b 2 >b 2c = - , (8) 

2 CL4 a 2 4 

and g = g c = b 2 ip c + b^ip 3 that is truncated by the liquid-gas coexistence 
boundary at the CEP b 2 = b 2e , if) = ip c (cf. case (a) in Figs. 8 and 9 of [10]). A 
detailed exposition of the Landau theory, with results for the phase boundaries 
and equilibrium values of <f> and if) in the various phases will be given in [18]. 



2 For B > 0, the mean-field correlation length of if) and hence the width of the 
interface region of classical kink solutions for if) are nonzero. The terms oc en and 
e 2 i are significant for relating the problem of critical adsorption of the a/3-phase at 
the a/3|7-interface to a wall problem. 
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To go beyond Landau theory, we use perturbation theory in combination with 
the RG.|3] Writing tp = tp Tei + tp, we expand TC about a reference value tp re{ , 
which we take as the mean-field value of tp at a reference point in the a(5 phase 
away from the critical line. This gives 



H[cP,tP}=H{0,tp rci ]+H'{<P,tP;tp 



ref 



(9) 



H'[(p,ip} = d d x 



Z fc=2,4 K 1=1 1 Z 



(10) 



with d 2 = a 2 +d2ifiei, a 4 = a 4 , h = -g+b 2 tp vei + b A (tp rei ) 3 , h = b 2 +3h(ip ie f) 
63 = 36 4 ipref, and 5 4 = 5 4 , where b 2 > 0. 
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If H! is taken into account by perturbation theory, the critical and first-order 
lines, and hence the CEP, get shifted. When studying the behavior near these 
lines, one must use their corresponding new locations that are compatible with 
the level of approximation. Suppose that a point (T c (g c ) , g c ) on the critical line 
is approached, which may be the CEP (g c = g e ). Let us ignore the tp 3 and 
tp 4 terms for the present. Then 7i becomes quadratic in tp, so that %p can be 
integrated out exactly. The resulting effective Hamiltonian 7i eT j[(p} is given by 
Hi[4>], with the replacements a 2 — > a 2 — (b\ d 21 /b 2 ) and Q 4 — > Cb^ — 

G^! /26 2 , PIUS 

a ^-independent term Jd d x /g(^2, &i) corresponding to the free energy of the 
Hamiltonian 7Yg[^] = Jd d x [(b 2 /2) ip 2 + bi tp}. As usual, we may presume that 
parameters such as d 2 , a 4 have a Taylor expansion in t = (T — T c (g c )) /T c (g c ) 
and 5g = (g — g c )/g c near (T c , g c ). Both H[0, ip re f\ as well as fc are regular at 
(T c , g c ). Hence the singular part of the total free energy results solely from H c s. 
That it has the usual scaling form near criticality can be demonstrated via a 
standard field-theoretic RG analysis, either of H e & directly, or of the equivalent 
model-C-type Hamiltonian from which it originated. Thus the critical behavior 
on the critical line is the same as at the CEP, provided our omission of the ip 3 
and tp 4 terms was justified. 

It is instructive to consider these nonlinearities first in the case d 2 \ = of a 
massive field ip decoupled from (p. Let <5tra[^] be the generator of transforma- 
tions of Hamiltonians T~C[tp] induced by the change of variable tp — > tp + 



0t»[*]W[$ = Jd d x 



5tp(x) 5tp(x)_ 



3 Implicit in our analysis is the well-founded assumption that the CEP, the critical 
line, and the first-order line g a will survive the inclusion of fluctuation corrections. 
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Choosing 



k bo 



^-i + 5 (0)^i^-3 



k = 3,4, 



(12) 



gives 

^ /cf^ = (25 tra [^ fc ] W G [^] +5 kA c 4 Jd d x, k = 3,4, (13) 



with c 4 = -36 4 [5(0)/26 2 ] 2 , where 5(0) = Jd d q/(2ir) d is a cutoff-dependent 
constant. Hence, at the Gaussian fixed point Hg[?P}, the -?/> 3 and ?/> 4 terms 
correspond to a redundant operator, and a redundant operator plus a con- 
stant, respectively. More generally, Wegner [22, Sec. III.G.2] has shown that 
any translationally invariant local operator can be represented as a redundant 
operator plus a constant at such a noncritical Gaussian fixed point. 

To generalize these considerations to the case d 2 i ^ with bi ^ 0, we insert 

CW^fe] W G [$ = <&tra[*fc] $ " / tfx f ^ 2 + 61) (14) 



into (13). The result tells us that, to first order in 6 3 , the ip 3 term is equiv- 
alent to shifts of a 2 , b 2 , and the constant part of 7i, plus a generated 2 ^/r 
contribution. Likewise, the ip A term corresponds (to order 64) to the gener- 
ation of 2 ?/> 3 and ?/> 3 contributions, and shifts of d 2 i and 61. Owing to the 
high naive dimensions 2(d — 1) and (5c? — 4)/2 of the operators 2 ?/> 2 and 
2 ?/r, we may trust that both produce only irrelevant corrections and hence 
may be dropped. Consequently, the effects of the -?/> 3 and -?/> 4 terms can be 
absorbed through shifts of the parameters d 2 , . . . , 61 of the Hamiltonian with 
&3 = &4 = 0, apart from irrelevant corrections. This means a change of the 
locations of the critical and first-order lines, and a corresponding adjustment 
of, e.g., the temperature scaling field. In summary, we arrive at an effective 
4 Hamiltonian H c s[<p], irrespective of whether the critical line or the CEP is 
approached. 

Next, we turn to the issue of the discontinuity eigenexponent y = d. Consider 
changes g c — > g c = g c + 5g, a 2c — > a 2c = a 2e + 5a 2 , away from the CEP 
such that the theory remains critical. As we have seen above, varying g (i.e., 
hi) alters the configuration-independent part of Ti, i.e., the coefficient, /x , 
of the 'volume operator' / d d x. Since fiQ trivially scales with the exponent d 
under RG transformations, it is formally relevant; but it does not contribute 
to the critical behavior at continuous phase transitions. Wegner [22] therefore 
proposed to call it 'special' scaling field. The above variation within the critical 
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manifold implies a change 5/j,q ~ (T c — T c ) of fiQ. This is the analog of the 
eigenperturbation with eigenexponent y = d found at the CEP fixed point in 
position-space RG calculations [5,6]. 

The above analysis of the symmetric CEP can be extended to the nonsymmet- 
ric case, in which the <fi — > — symmetry of H\4>, ip\ is broken [18]. As expected, 
the principal modifications are of a geometrical nature: the first-order surface 
bounded by the critical line is no longer confined to the h = plane, and the 
ordering field h and the temperature variable t mix in the scaling fields. 

Our analysis confirms that the singular part of the free energy of the aft and 
a + (3 phases has the usual scaling form anticipated in the phenomenologi- 
cal scaling theory, both on the critical line and at the CEP. Hence it is clear 
that the limiting forms (1) and (2) must hold and can be derived in a similar 
fashion as in phenomenological investigations [2,12]. We will therefore restrict 
ourselves here to a few remarks. The singularity (1) in ga{T) follows by ex- 
ploiting the equality of the free energies (grand potential) of the liquid (a[3, 
a, ft) and the spectator (7) phases at coexistence. To derive the behavior of 
the nonordering density p g = (■?/>) near the CEP, we have generalized the field- 
theoretic RG analysis [23] to the (^-dependent model-C-type Hamiltonian that 
result in the symmetric and nonsymmetric cases. In the nonsymmetric case ip 
does not only couple to the energy density (~ |^ 1 1 ~ Q: ) , but also to the order 
parameter, which produces the additional \t\^ singularity in (2). 
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